Opposite spin accumulations on the transverse edges by the confining potential 
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We show that the spin-orbit interaction induced by the boundary confining potential causes 
opposite spin accumulations on the transverse edges in a zonal two-dimensional electron gas in 
the presence of external longitudinal electric field. While the bias is reversed, the spin polarized 
direction is also reversed. The intensity of the spin accumulation is proportional to the bias voltage. 
In contrast to the bulk extrinsic and intrinsic spin Hall effects, the spin accumulation by the confining 
potential is almost unaffected by impurity and survives even in strong disorder. The result provides 
a new mechanism to explain the recent experimental data. 
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I. INTRODUCTION 



Recently, two experimental groups observed the trans- 
verse opposite spin accumulations near two edges of their 
devices in the presence of a longitudinal voltage biasji^ 
One experiment is on n-type GaAs's bar with a size about 
300 fim X 77/i77i,^ and the spin accumulation is detected 
by Kerr rotation spectroscopy. The other experiment is 
on a coplanar p-n junction light-emitting diode devicei^ 
Under a longitudinal bias, a circular polarization of the 
emitting light on two edges is detected. The directions 
of the polarization are opposite on the two edges, which 
suggests opposite spin accumulation on the two edges. 
Moreover, when the bias is reversed, the spin accumula- 
tions are reversed in the above two experiments. 

These experiments originally were motivated to mea- 
sure the predicted effects: the extrinsic and intrinsic spin 
Hall effects (SHE). The extrinsic SHE has been discov- 
ered about a few decades ago,^'* and it originates from 
the spin dependent scattering that deflects the spin-up 
and spin-down carriers towards the opposite edges of a 
sample. The intrinsic SHE is predicted first by Murakami 
et.al. and Sinova et.al in a Luttinger spin-orbit (SO) 
coupled 3D p-doped semiconductor^ and a Rashba SO 
coupled two-dimensional electron gas (2DEG)^, respec- 
tively. Recently, a number of sequent works have focused 
on this interesting effectiSiSiiS Nonetheless, the intrin- 
sic SHE still remains a controversy topic. The intrinsic 
spin Hall conductivity originally was pointed out to be 
universal in the clean bulk sample.- However some works 
have showed that the spin-Hall conductivity depends on 
the SO coupling strength and electron Fermi energy in 
generalii Moreover, it has been shown recently that the 
impurity plays an important role on the intrinsic SHE. In 
an infinite system, it has been shown that the spin Hall 
conductivity is very sensitive to disorder.* The effect van- 
ishes even in a very weak disorder limit when the vertex 
correction is considered. In a finite mcsoscopic ballis- 
tic system, the SHE can survive.^ The SHE and spin 
accumulations have been studied in the dirty or clean fi- 
nite mesoscopic samples by using the Landauer-Biittiker 
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FIG. 1: (Color online) Schematic diagram for the transverse 
effective potential Ve//(y) in the zonal 2DEG for the positive 

+Px- 



formalism and the tight-binding Hamiltonianii. These 
works show that the opposite spin accumulations indeed 
can generate on the transverse two edges in the finite sys- 
tem, and the SHE still presents below a critical disorder 
threshold.^ 

Although the spin accumulation observed in the two 
experiments appears to reflect the physics of the SHE, 
there is still a significant challenge to explain the ob- 
served experiment effect. The spin accumulation con- 
tributed from the extrinsic SHE has been shown to have 
the directions of spin polarization that are opposite to 
the experimental resultsi^ While the intrinsic SHE is pre- 
dicted to have the same directions of spin polarization as 
that observed in the experiments, it is expected to be 
scaled with length and eventually be vanished in Rashba 
spin orbit coupling systems in the presence of disorder. 
Since the experiments are not done in the ballistic region 
and the signal of spin accumulation has little dependence 
of the transversal length of the sample, it is not clear that 
how the intrinsic SHE explains all experimental results 
solely. 

In this paper, we give a new mechanism that can gen- 
erate the opposite spin accumulations in two edges under 
the longitudinal voltage bias. We first show the princi- 
ple of this new mechanism. Considering an infinite zonal 
2DEG with a confining potential on y-direction, which 

is described by the Hamiltonian H: H = ^^," -I- V{y), 
where Px/y is the momentum operator and m* is the ef- 



fective mass of electrons. V{y) is the confining poten- 
tial energy which is constant in the center and quickly 
increases as tending to the boundary. Based on the 
relativity effect, in the presence of the internal electric 
field E = VF(y)/e, there is a natural spin-orbit cou- 
pling term Vso = 2n%c^ cr • (E x p)j ^^i^^ Consider- 
ing that the corresponding electric field E = '^V{y)/e 
is perpendicular to the edge, only the element Ey of y- 
direction is non-zero and the spin-orbit coupling energy 
reduces into: Vso = - 2rJ'ici '^zPx-^V{y). So for the 
spin-down (cr^ =1 or —1) electrons, its effective potential 
Veffiy) = V{y) + Vsoiy) is lower than the one of spin- 
up electrons at the edge of y = for positive +px (see 
Fig.l). On the other hand, at the other edge of y = L, 
the effective potential Ve// of the spin-up electrons is 
lower for +px (see Fig.l). Thus, the spin accumulations 
on the two transverse edges form when electrons occupy 
the positive +px states under the positive longitudinal 
bias. When the bias is reversed, the electrons occupy the 
negative —p^ state and the spin accumulation reverses its 
sign. Therefore it produces the same spin accumulation 
as that observed in the experiments^^ In the present 
model, the spin accumulation originates completely from 
the structure confining potential, so it is not affected by 
the impurity and the dcphase. In a word, the structure 
confining potential can also induce the opposite spin ac- 
cumulation, which can be the origin of the experimentally 
observed spin accumulation. 

The paper is organized as follows: in section II, we will 
mention and solve the model in detail. The results and 
discussions are in section III. Finally, a brief summary is 
given in section IV. 



II. THE MODEL AND SOLUTION 



The Hamiltonian of the zonal 2DEG can be written as: 



where 7 = .,^,^^"2 ■ 



In the experiment, the length L along 
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a^Px-rViy). (1) 



Here the first term is kinetic energy, the second term 
is the potential energy, and the third term is from the 
spin-orbit coupling energy due to the boundary confining 
potential V{y) as mentioned in the introduction. 

Due to the fact [px,H] = and [az,H] — 0, k^ and 
az are the good quantum numbers, the eigenstates in 
such a structure can be written in the form \I'(x,y) = 
4>nk^ {y)exp(ikxx) with the dispersion relation E — enk^ + 
h^k'^/2m*. The index n numbers the different subbands 
with the wave-function (l)nk^ (y) in the y-direction and 
a energy enk^ , which are functions of the wave vector 
kx- The transverse wave- function 4>nk^{y) satisfies the 
equation: 



2m* ay 



(f'nk^iy) = f^nk^4>nk^{y), 



the transverse y-direction is very long [e.g. in the order 
of tens of micron in Ref.(l)] while the confining potential 
is only limited in several atom's layers. So we model this 
potential V{y) as a square potential well, i.e. V{y) = 
for < y < L and V{y) = V for others y. In this case, 
Vso = ~l'^z-^V{y) becomes the 5-function. The Hamil- 
tonian can be solved analytically, and the Schrodinger 
equation Eq.|21) for the spin-up electron reduces into: 



^2m* 



V{y) + 7( S{y) - S{y - L) ) ]0(y) = e0(y). (3) 



The spin-down electronic Schrodinger equation is same 
to Eq.Q except 7 — > —7. The wave function 4>(y) in 
Eq.(|3J) can be written as: 



^{y) 



Aef^y, y < 

sin{kyy + 0), < y < L 
Be-Py, y> L 



where ky = \/2m*e/h, j3 = ^j2m*{V — ()/h, and A, 
B, and 6 are the constants to be determined by the 
boundary conditions. Here the boundary conditions 
are 0(y)|y=o-/L- = '/'(y)ly=o+/L+ and (/>'(y)|y=o+/L+ - 
(t)'iy)\y=o-/L- = ±i'^m*-//h'^)(j){y)\y^o/L, which lead to 
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sin{kyL + 0) =Be-'^^ 

kyCosB - A(3 ^ 2m*-/ A 

PBe-^^ + kyCos{kyL + 6) = 2m*-fBe- 



■fSL 



Solving the equation set, we get e„fc^ and 4>nk^{y)- Se- 
quentially the spin-up electron probability distribution 
Pnk^,l{y) = |0nfc^(2/)P is obtained. The same calcula- 
tion can be done for the spin-down electron. Due to the 
system is universal along the x-direction, and the spin 
accumulation is independent of x, Pnk^,iiiy) completely 
describes spin distribution. 

Although the square confining potential V{y) model 
is solvable analytically, the potential in the real system 
is not abrupt. An gradual change from bottom to top 
close to the interface is expected. For this reason, we 
consider the real parabolic confining potential V{y) and 
solve the system numerically by using the tight-binding 
Hamiltonianii. In addition, we also study the disorder 
effect on the spin accumulations. In the tight-binding ap- 
proximation, the Hamiltonian in Eq. ^ , which is related 
to y-direction, can be written as the following discrete 
lattice version: 



H 



i,a 



{e^ -h K: + o-zVso,i)al^aicr + V^ ^ ^L^jcr, (4) 
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where a =t,i (or ±1) is the spin index in z-direction, 
and t = h?/2m*a'^ represents the hopping matrix ele- 
ment with the lattice constant a. The confining poten- 
tial V{y) is assumed to be parabolic: V{y) — V 250° 



(2) for a < y < 5a, V{y) = for 5a < y < L — 5a, and 



V{y) = Vi^^±|^ ior L - 5a < y < L. For a clean 
system, the onsite energy e^ = 0, and e^ is randomly dis- 
tributed between [—W/2, W/2] for the dirty system. The 
Hamiltonian in Eq. I^J can easily be solved by numerically 
calculating the eigenvalues and eigenstates of the Hamil- 
tonian matrix. Due to decoupling between the different 
spin states in the Hamiltonian, we can solve the eigen- 
values and eigen-wavefunctions separately for the spin- 
up and spin-down electrons. Whereafter, the the spin-up 
and spin-down electron probability distribution Pnk^cdi^) 
in the subband n and the longitudinal momentum hkx is 
obtained straightforwardly. 

While the longitudinal bias is zero, the spin accumu- 
lations S{y) is zero everywhere, because of the existence 
of the time-reversal invariance. On the other hand, when 
a bias Vuas is added, the spin accumulations S(y) will 
emerge. Consider the device under the positive bias Hms 
and at the zero temperature, the -\-kx states with its 
energy between Ef — Vuas/'i, and Ef + Vuas/'i, are oc- 
cupied by electrons, while the states for the negative 
— fcj. are empty, here Ef is the Fermi energy. Then, 
under the small bias and taking the linear approxima- 
tion, the spatial density distribution P|x(i) of the spin- 
up (down) electrons along y-direction for the unit bias 
is P^iH) = J2nP'^(^kjPnk^,'\ii), where Pn{EkJ is the 
density of state in the subbands n with Ek^ = Ef — Cnk^ , 
and the sum is over all subbands n with its cut-off en- 
ergy lower than Ef. The spin accumulation density 
and charge density in the linear bias can be obtained 



linn/„„,^c 



S{ia) 



|(Pt(0 - n(*)) and 



as: Ps{i) = iimv^i^^^o y^ 

PeW-e(PT(z) + PxW). 

In the numerical calculation, we choose the realistic 
parameters as ones in the experiment: ^^ the electronic 
effective mass m* = 0.05me and the Fermi energy Ef = 
O.ley, in which the corresponding electron concentration 
is approximately n2D — 10^^ cm^^. The energy unit is 
set to leV, and the length unit is Inm for the analytic 
model, or la.u. for the tight-banding model. The lattice 
constant o in the tight-banding model is around 0.196n?7i. 



III. NUMERICAL RESULTS AND DISCUSSION 

First, we study spin accumulation in the clean system. 
In the present system, the spin accumulation density 
Psiy) depends on the transverse position y, and it is in- 
dependent of the longitudinal position x. Fig. 2 and Fig. 3 
show the spin accumulation density Ps{y) and charge 
density Pe (y) versus y for the case of the square and the 
parabolic confining potential V{y), respectively. First 
of all, the opposite spin accumulations at z-direction in- 
deed are generated near the two edges regardless of the 
square or the parabolic potential. If the bias is reversed, 
the electrons occupy the negative —k^ states instead of 
the positive +kx states and the spin accumulations also 
are reversed. These results are consistent with the ex- 
perimental results;!^ The opposite spin accumulations 
here obviously originates from the confining potential as 
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FIG. 2; (Color online) The transverse distribution of the spin 
and charge density, Ps{y) and Pe{y), vs y for the square 
confining potential. Left panel: L = 29.3nm, right panel: 
L — 58.6nm. 
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FIG. 3: (Color online) The transverse distribution of the spin 
and charge density, Psiy) and Peijj), vs y for the parabolic 
confining potential. Left panel: L = 150a, right panel: L = 
300a. 



mentioned in the introduction because there is no other 
interaction except for the potential V{y). Second, the 
spin accumulations mainly are near the two edges, and 
it is small and has oscillation in the bulk. The oscilla- 
tion is expected due to the existence of Fermi surface. 
The period is given by 27r/fci?. For a fixed Fermi en- 
ergy (e.g. Ef — Q.leV), the wider the width L is, the 
more the subbands below Ef will be. At the mean time, 
the oscillation times of Ps and Pe are more and the os- 
cillation amplitude are smaller. So the bulk Ps almost 
vanishes and Pe approaches constant at large L. But 
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FIG. 4: (Color online) The spin polarization Ps/Pe vs the 
transverse site y for the square confining potential. Panel 
(a),(b),(c),(d) correspond to L = 14.6nm, L = 29.3nm, L = 
58.6nm, and L = 117. 3nm, respectively. The insets show the 
detailed distribution near the edge of j/ = in the central 
panel. 



FIG. 5: (Color online) The spin polarization Pa/Pe vs the 
transverse site y for the parabolic confining potential. Panel 
(a),(b),(c),(d) correspond to L = 75a, L = 150a, L = 300a, 
and L = 600a, respectively. The insets show the detailed 
distribution near the edge of y = in the central panel. 



the spin accumulations near the edge, including the in- 
tensity and the location, is almost independent with L. 
Third, we discuss the characters of the spin accumula- 
tions Ps (y) as a function of the strength V of the confin- 
ing potential V{y). The characters are slightly different 
for the square and paraboHc potentials. While V — 0, 
Ps{y) — for both the square and parabolic potentials. 
With V increasing, Psiy) emerges. For the square po- 
tential, Psiy) quickly arises in the beginning. Around 
V = 5Ef, Ps{y) has reached saturated value. Thereafter 
the value almost does not change with further increas- 
ing V (see Fig. 2a, b). For the parabolic potential, Ps{y) 
increases slowly. Around V — 8eV, it saturates. For 
comparison, we also show the charge density Pe{y) in 
Fig. 2 and Fig. 3. Here Pe{y) > in any position y, and 
Pe{y) is symmetric while Psiy) is asymmetric. 

Next, we discuss the spin polarization Ps{y)/Pe{y)- In 
the central panel of Fig. 4 and Fig. 5, we plot the trans- 
verse distribution of the spin polarization for the square 
and parabolic confining potential V{y), respectively. It is 
clearly shown that the opposite spin polarization emerges 
on the transverse two edges whatever L is set. For ex- 
ample, in the Fig. 4 and 5, the transverse width L is set 
to 14.6n?7i, 29.Snm, 58.6nm and 117.3nm which is very 
different, but the spin polarization distributions near the 
edge are almost identical. In addition, the spin polariza- 
tion on the transverse edge is close to a linear function 
as V while Ps{y) is hardly affected by the strength of V 
(see Fig. 2 and 3). The detailed distributions of the spin 
polarization near the edge oi y — are magnified and 
shown in the insets (in Fig. 4 and 5). The distribution 
range for the parabolic potential is slightly wider than 
that for the square potential because the variation range 
of the parabolic potential V{y) is wider than that of the 
square potential V{y). 



The above result is obtained in a clean system. 
In the following, we study the spin accumulation 
in the dirty system. Fig. 6 displays the transverse 
distribution of Psiy), Pe{y) and spin polarization 
Ps{y)/Pe{y) with different disorder strength W = 
0.lEf,lEf,2Ef.5Ef,10Ef. In these calculation, P,(y) 
and Pe (y) are obtained by averaging over up to 5000 re- 
alizations of disorder. From Fig. 6, we can see that the 
spin accumulation Psiy) and the spin polarization near 
the two edges are almost unaffected by the disorder, even 
the disorder W is very strong (e.g. W — WEf). In fact, 
the spin accumulation in the present device originates 
from the confining potential near the edge, where the ef- 
fective potentials for spin-up and spin-down electrons are 
different (see Fig.l). Intuitively, the spin polarization is 
expected to be unaffected by the disorder as well as the 
dephase. Additionally, with the increasing of the disorder 
W, the amplitudes of the oscillation of Pe{y) and Ps{y) 
in the bulk are slightly reduced (see Fig. 6a, b), and their 
fluctuations are increased linearly. 

To compare our numerical results with the 
experimenlii, we calculate the value for spin polar- 
ization and accumulation. From the experimental 
data, the spin density at the peak can be estimated: 



Ps « (1.5- -4.2) X IQ- 



eV , and the spin polar- 



ization is (1.0 2.8) X 10 ** for a thickness h = 0.9fj,m. 

From figures (e.g. fig. 2 and fig. 4), our calculation show 



that the peak spin density Pg ~ 0.3 x 10 



'eV 



and the spin polarization is 0.8 x 10^''. The value of the 
spin polarization is comparable with the experiment. 
Due to the big thickness in the cxpriment, the spin 
accumulation Ps in our calculation seems to be several 
times smaller than that of the experiment. Indeed, a 
more precise quantitative calculation perhaps requires 
to consider the spin orbit coupling in the bulk and spin 
relaxation effect at the boundary. However, our simple 
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FIG. 6: (Color online) The charge density distribution P^ 
(a), the spin density distribution Ps (b), and spin polarization 
Ps/Pe (c) vs the transverse site y for the different disorder W, 
respectively. The other parameters are: i?/=0.1eV, L = 150a, 
and V = 8eV. 



mensional electron system. Due to the strong structure 
confining potential in the boundary, the induced spin- 
orbit interaction leads to the opposite spin accumulation 
on the two transverse edges under the longitudinal volt- 
age bias. The spin polarized direction can be reversed 
while the bias is reversed. The intensity of the polar- 
ization is also proportional to the external longitudinal 
voltage bias. These results are consistent with the re- 
cent experiment. Moreover, the experimental test of the 
new mechanism can be easily performed in future exper- 
iments. Unlike the extrinsic and intrinsic SHE, the spin 
accumulations in the present mechanism are hardly af- 
fected by the disorder and dephase, and can exist even 
in the strong disorder system. 



model indeed produces the same order of magnitude as 
one measured in the experiments. 
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